/ 

4 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  This  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  The  lime  for  reviewing  instructions,  searching  existing  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  esiimale  or  any  other  aspect  of  this  collection  ol 
inlormation,  including  suggestions  for  reducing  the  buiden,  to  the  Department  of  Defense,  Executive  Services-  and  Communications  Directorate  (0704-0 IflSl  Respondents  should  ba  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penally  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  QMB 
control  number, 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  IDD-MM-YYYY)  2.  REPORT  TYPE 

22-10-2007  Journal  Article 

3.  DATES  COVERED  (From  ■  To! 

4.  TITLE  AND  SUBTITLE 

Why  Does  the  Secchi  Disk  Disappear?  An  Imaging  Prospective 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

060243 5N 

6.  AUTHORISI 

Weilin  Hou,  Zhongping  Lee,  Alan  D.  Weidemaan 

Sd.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

73-6867-07-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stcnnis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/JA/7330-06-7053 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  ANO  ADDRESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  22217-5660 

10.  SPONSOR/MONITOR  S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR'S  REPORT 

NUMBERIS) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  widely-used  Secchi  disk  method  is  re-examined  from  the  modulation  transfer  aspect  Namely,  by  assuming  a  volume  scattering  function  and 
applying  small  angle  scattering  approximation,  we  show  that  the  Secchi  depth  and  horizontal  visibility  can  be  determined  using  the  water  modulation 
transfer  function  and  the  corresponding  spatial  frequencies  associated  with  the  disk.  A  basic  equation  of  Secchi  disk  is  reached  that  is  comparable  to 
the  radiative  transfer  approach,  in  that  the  Secchi  disk  is  reached  that  is  inversely  proportional  to  the  attenuation  coefficient  (c).  With  typical  values  for 
parameters  applied,  we  demonstrate  that  the  modulation  transfer  technique  produces  a  horizontal  visibility  range  of  about  4,8/c,  which  is  inline  with 
previous  studies  The  improvement  lies  in  the  fact  that  the  current  approach  correctly  addresses  the  repsponse  of  all  spatial  frequencies  according  to 
the  modeled  optical  transfer  function  of  the  water  In  terms  of  Secchi  disk  theory,  the  current  approach  helps  to  understand  the  effect  of  disk  size  as 
well  as  the  role  of  scattering  of  the  Secchi  disk  depth.  The  approach  presented  provides  an  understanding  of  Secchi  disk  disappearance  by  showing 
that  the  disk  is  moved  away  from  the  ovserver,  the  spatial  frequencies  corresponding  to  the  disk  size  increase,  while  the  modulation  transfer  dampens 
contrast  at  an  increased  rate 

15.  SUBJECT  TERMS 

atmospheric  and  oceanic  optics,  imaging  systems,  vision,  color  and  visual  optics 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

me 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Weilin  Hou 

a.  REPORT 

b,  ABSTRACT 

c,  THIS  PAGE 

Ur 

PAGES 

12 

Unclassified 

Unclassified 

Unclassified 

UL 

19b.  TELEPHONE  NUMBER  (Include  area  code 1 

228-688-5257 

Standard  Form  298  {Rev.  8/98) 

Proscribed  by  ANSI  Std. 


Why  does  the  Secchi  disk  disappear?  An  imaging 

perspective 

Weilin  Hou,  Zhongping  Lee,  Alan  D.  Weidemann 

Naval  Research  Lalxtratory.  Code  7333 ,  Stennis  Space  Center ;  MS  39529 
hou  (^nrlsst  \  navy,  mil 


Abstract:  The  widely-used  Secchi  disk  method  is  re-examined  from  the 
modulation  transfer  aspect.  Namely*  by  assuming  a  volume  scattering 
function  and  applying  small  angle  scattering  approximation,  we  show  that 
the  Secchi  depth  and  horizontal  visibility  can  be  determined  using  the  water 
modulation  transfer  function  and  the  corresponding  spatial  frequencies 
associated  with  the  disk.  A  basic  equation  of  Secchi  disk  is  reached  that  is 
comparable  to  the  radiative  transfer  approach,  in  that  the  Secchi  depth  is 
inversely  proportional  to  the  attenuation  coefficient  (c).  With  typical  values 
for  parameters  applied,  we  demonstrate  that  the  modulation  transfer 
technique  produces  a  horizontal  visibility  range  of  about  4.8/t\  which  is 
inline  with  previous  studies.  The  improvement  lies  in  the  fact  that  the 
eunrem  approach  correctly  addresses  the  response  of  all  spatial  frequencies 
according  to  the  modeled  optical  transfer  function  of  the  water.  In  terms  of 
Secchi  disk  theory,  the  current  approach  helps  to  understand  the  effect  of 
disk  size  as  well  as  the  role  of  scattering  on  the  Secchi  disk  depth.  The 
approach  presented  provides  an  understanding  of  Secchi  disk  disappearance 
by  showing  that  as  the  disk  is  moved  away  from  the  observer,  the  spatial 
frequencies  corresponding  to  the  disk  size  increase,  while  the  modulation 
transfer  dampens  contrast  at  an  increased  rate. 
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1,  Introduction 

Ocean  optics  can  trace  its  origin  to  understanding  and  determining  underwater  visibility.  The 
simple,  widely-used  Secchi  disk  method,  dating  back  to  1865,  Is  still  used  by  biologists  and 
physicists  as  a  quick  measurement  of  water  clarity  in  natural  waters*  Its  deployment  is  simple: 
one  lowers  the  traditionally  white*  circular  disk  about  30cm  in  diameter,  from  above  the  water 
into  the  water  column,  and  determines  the  point  at  which  it  disappears  from  sight.  The  depth 
at  which  the  disk  can  no  longer  be  discerned  against  the  background  is  defined  as  the  Secchi 
disk  depth  (Z5D).  Measurement  protocols  include  measurement  of  the  depth  at  which  the  disk 
disappears,  and  then  re-appears  and  then  using  the  average  of  both  values  to  eliminate  emirs 
imposed  by  possible  residual  memory  effects,  and  by  different  operators.  Secchi  depth  is  often 
used  as  a  convenient  way  of  quantifying  water  clarity.  It  could  also  be  used  horizontally,  to 
help  measure  diver  visibilities,  or  to  avoid  issues  associated  with  its  vertical  deployment  (ie 
sunglim,  surface  waves,  bottom  patchiness  etc),  especially  in  very  shallow  areas. 
Preisendorfcr  [1,2]  thoroughly  summarized  and  reviewed  the  science  behind  the  Secchi  disk 
approach  from  the  radiative  transfer  perspective,  and  warned  against  potential  over-usages  of 
its  application  in  what  he  termed  “Secchi  Disk  madness.”  Even  though  new,  more  rigorous 
electro-optical  instruments  are  available  to  measure  optical  properties,  the  Secchi  depth  still 
has  proven  to  he  an  inexpensive,  dependable,  and  efficient  measure  of  the  water  darily  that 
has  been  related  to  diver  visibility,  water  quality,  and  biological  activities  [2,  3].  The  vast 
database  of  available  Secchi  depth  accumulated  over  the  years  throughout  the  entire  world, 
nonetheless,  adds  the  incentive  to  better  understand  the  science  behind  this  simple  approach. 
The  implications  of  the  Secchi  disk  method  on  underwater  imaging  issues  such  as  prediction 
of  diver  visibility  for  mine  identification  also  prompt  further  study  beyond  the  traditional 
radiative  transfer  approach. 

For  most  theoretical  treatments  of  Secchi  disk  visibility,  the  common  definition  of  the 
visibility  contrast  (Weber  contrast)  is  used  [1,2,  4].  We  shall  follow  conventional  notations 
cited  and  briefly  outline  the  theory.  Details  can  be  found  in  above  mentioned  references.  The 
visibility  contrast  Cv  is  given  as  the  normalized  difference  in  brightness  (radiance, 

L(#,0,  z) )  between  the  target  (Secchi  disk) 


c 


(i) 


where  subscripts  T  and  B  denote  the  target  (disk)  and  the  background  respectively.  As  the 
viewing  range  decreases  to  zero  (without  the  medium  attenuation  and  scattering),  the  above 
gives  the  inherent  contrast  (CTO).  For  horizontal  orientation  and  viewing,  as  the  disk  is  moved 
further  away  (increasing  z)  it  can  be  shown  that  the  apparent  contrast  (C vz)  decreases 
exponentially  as  function  of  the  medium’s  attenuation  [2,  5]: 


Cyz  —  CVf)e 


(2) 


where  c  is  the  photopic  beam  attenuation  coefficient,  which  is  the  attenuation  of  the  natural 
light  spectrum  convolved  with  the  spectral  responsivity  of  the  human  eye  (photopic  response 
function).  For  most  viewing  conditions  this  can  be  approximated  by  a  monochromatic  beam 
attenuation  near  the  peak  of  the  human  eye  sensitivity  such  as  at  532nm  [4].  When  the  disk  is 
lowered  vertically  into  the  water  column,  one  must  take  into  account  that  the  illumination  at 
the  disk  and  background  is  being  influenced  by  the  surface  radiance  and  attenuation  of  light, 
therefore  the  above  equation  is  modified  by  the  diffuse  attenuation,  and  becomes  [2, 5] 
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Cvz=Cvoe-(c+K)',  (?) 

where  K  is  the  downward  diffuse  attenuation  coefficient  (ignoring  the  viewing  angle 
dependence  for  simplicity).  The  above  (Eq,  (3))  is  what  Preisendorfer  referred  to  as  the  basic 
equation  of  Secchi  disk  science-  Or  in  a  different  form  when  the  contrast  of  the  target  and 
background  converge 


c  +  K=- 


1n(C,) 


CL  = 


J SD 


' V7. 


"VO 


(4) 


where  Q.  is  the  limiting  contrast  at  disappearing  disk  depth  (ZSD)  [4],  CL  has  been  shown  to  be 
a  function  of  the  disk  size,  but  more  importantly  the  adaptation  luminance  at  the  disk  location 
[6].  The  contrast  threshold  has  been  found  to  vary  from  0.02  or  higher  with  low  adaptation 
luminance,  to  0.008  [7]  or  even  0.002  at  high  luminance  [6].  Unless  mentioned,  all  parameters 
are  associated  with  photopic  quantities  hereon.  It  can  be  shown  that  such  simplification  does 
not  invalidate  derivations  that  will  follow  [4].  The  spectral  and  angular  dependence  of 
radiance  distribution,  as  well  as  the  effects  of  surface  glint,  are  also  neglected  tor  simplicity. 

When  a  black  Secchi  disk  is  used  instead  of  the  traditional  white  one,  studies  suggest  that 
the  same  theory  holds  with  less  variability  [7],  and  therefore  the  black  disk  can  be  used  as  a 
robust  measure  of  underwater  visibility  [4].  Other  modifications  of  the  disk,  such  as  one  with 
alternating  hlaek  and  white  quadrants,  are  also  used  by  many  researchers,  especially  those 
monitoring  lakes.  Results  show  visibility  ranges  measured  with  a  white/black  quadrant  disk 
are  similar  to  those  obtained  with  the  single  color  black  disk  [8].  This  has  important 
implications  from  the  imaging  perspective  in  that  the  white/black  quadrant  disk  can  be  viewed 
as  a  pseudo-resolution  chart,  with  a  separation  distance  about  half  the  size  of  the  disk 
diameter.  This  makes  sense  if  we  consider  a  disk  with  half  of  its  side  painted  black  and  half 
painted  white,  which  can  be  approximated  to  a  bar  pattern  with  a  cycle  of  d .  Splitting  the  area 
again  to  the  quadrant  Secchi  disk  reduces  the  size  of  the  pattern  by  half  to  approximately  a  df 2 
cycle. 

Intuitively,  since  the  disappearance  of  Secchi  disk  is  observed  by  a  human  operator,  it  is 
fitting  to  examine  the  issue  from  the  imaging  perspective.  In  this  regard,  underwater  visibility 
can  be  thought  of  as  the  resolution  reduction  of  specific  spatial  frequencies  over  distance. 
This  is  usually  studied  by  examining  a  standard  resolution  chart  such  as  USAF  1951  and 
known  environmental  transmittance  characteristics.  To  understand  how  the  visibility  contrast 
fits  the  image  theory,  it  is  prudent  to  review  the  background  briefly.  Details  can  be  found  in 
many  image  processing  books  such  as  [9]. 

Generally  speaking,  a  2nii mention al  image  of  an  object  is  basically  the  combination  of 
original  signal,  fi[x,y).  convolved  by  the  imaging  system  response  h(xTyl  integrated  over 
sensor  space  E: 

g(x,y)  =  j\f{xi,yi)h(x-xi,y-yi)dxldyl ,  (5) 


or, 

g(x,y)  =  f(x,y)®h(x,y),  (6) 

here  ®  denotes  2-D  convolution,  and  h(x,y)  is  the  system  response  to  a  point  source,  or  the 
pin nt -spread  function  (PSF),  The  system  response  includes  those  from  hoth  the  imaging 
system  itself  as  well  as  the  effects  of  the  medium  (water  in  our  case). 

Mathematically,  it  is  easier  to  manipulate  the  above  relationship  in  the  frequency  domain 
as  the  convolution  operator  becomes  simple  multiplication.  Applying  Fourier  transform,  the 
above  relationship  becomes 
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G(u%  v)  =  F{u%  v)H(u t  v) ,  (7) 

where  u,  v  are  spatial  frequencies  and  G*  F,  H  are  Fourier  transforms  of  g *  /  and  h 
respectively  .  The  Fourier  transfer  of  /i*  for  example*  is  in  the  following  form: 

ff(«,v)  =  J  \h(x,y)e-ilt<UI^)dxdy.  (8) 


The  system  response  function  H *  also  referred  to  as  the  optical  transfer  function  (OTF)*  is 
the  Fourier  transform  of  the  PSF.  The  magnitude  of  the  OTF  is  the  modulation  transfer 
function  (MTF).  The  MTF  describes  the  contrast  response  of  a  system  at  different  spatial 
frequencies*  and  when  the  phase  information  is  of  little  concern,  it  is  a  sufficient  measure  of 
the  power  transfer.  By  definition*  therefore*  the  MTF  can  he  measured  by  the  contrast  of 
sinusoid  or  bar  patterns  of  corresponding  spatial  frequencies  [9*  10|*  and  is  the  focus  of  this 
study. 

Notice  that  the  above  MTF  term  H(u>v )  is  the  total  system  response.  Therefore  if  one 
views  the  complete  path  from  target  to  the  bottom  of  eyes  or  the  recording  CCD  plane*  in 
many  cases,  the  MTF  can  he  the  effect  of  multiple  individual  components.  In  the  frequency 
domain,  the  MTF  can  be  expressed  by  the  direct  product  of  each  component,  for  instance,  the 
optical  system  itself,  and  the  medium  (plus  any  other  factors  when  applicable): 

H(u,v)  =  (9) 


The  system  response  H3y5trm(utv)  can  be  p re-determined  and  calibrated  to  remove  any 
significant  errors*  and  in  most  cases,  does  not  vary  with  imaging  conditions.  The  above 
formulation,  which  emphasizes  the  validity  of  the  separation  of  the  system  and  the  medium*  is 
significant  in  our  analysis.  Furthermore*  one  should  pay  special  attention  to  the  band-limiting 
characteristics  imposed  by  such  as  a  camera  system's  field-of-view,  and  Nyquist 

sampling  frequency  limits  imposed  by  the  CCD  resolution  [9],  For  a  known  object*  by 
examining  the  characteristics  of  the  medium  evident  in  //Wi/im7I(w*v),  one  could  theoretically 
predict  the  exact  outcome  image  within  that  environment  by  responses  from  all  of  the 
individual  spatial  frequencies.  In  contrast*  one  can  inversely  derive  the  detailed  information 
of  the  object  from  the  outcome  image  via  inverse  Fourier  transform,  which  is  the  goal  of 
underwater  mine  detection  and  other  applications  including  target  recognition  and  tracking. 
Needless  to  say,  the  presence  of  various  noises  (such  as  scattering  or  surface  fluctuations) 
complicates  these  through-the-sensor  techniques  and  is  beyond  the  scope  of  this  paper. 
Instead*  the  focus  here  is  on  the  optical  properties  which  influence  the  modulation  transfer 
functions  for  the  underwater  environments  as  they  relate  to  the  disappearance  of  the  Secchi 
disk.  To  better  illustrate  the  main  problem  from  the  imaging  point  of  view*  we  choose  optimal 
viewing  conditions  associated  with  Secchi  depth  measurement:  minimal  influence  from 
surface  waves,  sunglinh  shadows,  and  bottom  reflectance.  For  convenience  we  focus  on  the 
horizontal  visibility  range  with  ZSD  used  to  generally  refer  to  the  disappearance  range.  While 
we  focus  on  the  horizontal  case*  the  vertical  depth  can  be  derived  similarly  from  the 
horizontal  results  by  simply  adding  photopic  diffuse  attenuation  K  to  the  c  values  Eqs.  (2X3). 

2,  Secchi  theory  based  on  modulation  transfer 

Recall  that  the  modulation  transfer  function  measures  the  modulation  or  contrast  of  a  pattern 
after  transmission  through  the  system.  The  modulation  can  be  obtained  by  measuring 
maximum  and  minimum  values  of  the  sinusoid  (or  bar)  pattern  of  spatial  frequency  ip  [9*  10], 

*  (l0) 
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where  and  Smm  denote  the  maximum  and  minimum  ‘‘brightness1”  values  of  the  pattern  at 
each  spatial  frequency  ^  at  range  Z.  These  terms  can  be  related  to  the  luminance  level  of  the 
image,  or  the  radiance  detected  by  the  image  sensor. 

By  normalizing  to  the  input  or  initial  modulation  (assuming  uniform),  which  is 


M(p,0)  = 


^(^0)-^m,n(^0)  =M 

SMW)  +  Smi¥,0)  °* 


(ID 


the  modulation  transfer  function  H{ \p,Z)  can  be  expressed  as 


H(y.z)  = 


M0 


(12) 


When  considering  black  versus  while  patterns  the  source  modulation  is  unity,  M0  (=1)  can 
be  omitted.  Operationally,  normalization  can  be  applied  by  dividing  the  low  frequency  value 
of  the  system  M0  before  the  transfer*  Note  that  while  the  above  contrast  definition  (also  known 
as  Michelson  contrast)  is  different  from  the  traditional  contrast  equation  used  in  Secchi  disk 
theory,  when  the  Secchi  disk  disappears  the  information  content  is  the  same  in  both 
formulations. 

For  circular  symmetrical  response  systems,  such  as  the  isotropic  volume  scattering  type 
found  in  the  seawater,  the  corresponding  2-dimensional  transforms  found  in  Eq,  (8)  can  be 
reduced  to  a  one-dimensional  Hankel  (Fourier-Bessel)  integral,  such  as  for  H: 

H(y/yz)  =  2x  ^  JQ{2n&y/)h(0yz)6ti0 .  (13) 

8=Q 


Wells  (II]  applied  small  angle  approximations  to  the  above  and  derived  a  robust 
underwater  modulation  transfer  model  which  is  briefly  outlined  below*  By  separating  the 
exponential  decay  effect  with  distance  due  to  the  medium,  the  MTF  of  the  medium  in  Eq.  (13) 
can  be  expressed  as 

H(i//,z)  =  e'0i¥)Z ,  (14) 

where  D{tp)  is  the  decay  transfer  function  (DTF)  and  is  independent  of  the  range  of  detection. 
This  provides  a  method  to  compare  measurements  at  different  ranges  for  consistency. 

By  using  a  thin  slab  model  with  the  small  angle  scattering  approximation,  the  decay 
function  can  be  written  as  the  photopic  beam  attenuation  c,  less  the  light  scattered  back  into 
the  acceptance  cone,  which  is 

D(v)  =  c-S(v)>  (15) 


where 


Sty)  =  2x  jjo(2xOr)0(0)6d0.  (16) 

and  fl(ft)  is  the  volume  scattering  function.  Total  scattering  coefficient  h  is  obtained  via 

b  -  2k^  f3{0)  sin  &J0  *  (17) 

In  an  effort  to  derive  a  closed-form  solution,  Wells  [111  assumed  a  scattering  function  with 
the  following  analytical  form. 
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(18) 


P{6)~  27T{dt+d2)m' 

where  00  is  related  to  the  mean  square  angle  (MSA).  When  compared  with  coastal 
measurements  by  Petzold  [12]  (Fig,  I  )*  this  function  approximates  the  behavior  of  the 
scattering  in  small  angles  reasonably  well  in  most  cases,.  More  discussion  about  these  curves 
will  follow  in  the  next  section.  Combining  Eqs,  (I5MI8),  the  DTP  of  the  seawater  can  be 
expressed  as 

D(y)  =  c-S(y/) 


bQ-e -***) 


(19) 


curves  are  normalized  at  T  for  comparison. 

Equations  (14)  and  (19)  above  reveal  the  relationship  between  the  water  MTF,  the  total 
attenuation,  the  spatial  frequencies  involved  (which  can  be  related  to  the  physical  dimensions 
of  the  target  or  Secchi  disk),  and  the  depth  or  horizontal  range  at  which  the  Secchi  disk 
disappears.  For  example,  assume  single  scattering  albedo  a>o=  6/c= 0.8,  00=0.03,  cZ= 4,8,  the 
MTF  of  the  water  plotted  in  Fig,  2  shows  a  large  and  rapid  reduction  with  increasing  angular 
frequencies,  following  Eq.  (12).  It  is  clearly  shown  that  finer  details  (higher  frequencies)  are 

lost  first  Notice  that  the  term  only  affects  very  low  spatial  frequencies. 

From  a  qualitative  MTF  perspective  it  is  clear  that  the  disappearance  of  the  Secchi  disk  as 
it  is  moved  away  from  the  observer  is  the  result  of  reduced  resolution  by  the  system  response 
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function,  or  MTF,  at  the  spatial  frequency  related  to  the  disk  size  and  range.  Here  for 
simplicity,  we  exclude  the  effect  of  imaging  system  itself  (the  human  eye  in  this  situation),  as 
well  as  that  of  the  air-sea  interface  when  considering  vertical  application,  and  consider  only 
the  effect  of  water  transmission. 


Fig.  2.  MTF  of  water  at  different  mean  square  scattering  settings  Eq,  (19),  with  optical  length 

DZ=4.8.  The  curve  with  circles  demonstrate  the  MTF  without  €  term  in  Eq.  (19).  Sec 

text  fur  details. 

When  the  limiting  contrast  threshold  CL  is  reached,  the  Secchi  disk  can  no  longer  be  seen 
by  the  human  operator  The  resulting  effective  spatial  frequency  (Tsn)  reflects  all  frequencies 
related  to  the  disappearance  of  the  disk.  These  include  those  related  to  the  actual  disk  size 
(low  frequencies)  extending  to  those  associated  with  the  disk  edge  (high  frequencies).  Recall 
our  earlier  discussion  regarding  the  alternating  black  and  white  Secchi  disk  as  a  pseudo 
resolution  chart  with  a  separation  distance  of  dll ,  the  first  order  approach  to  approximate  ¥sd 
can  he  expressed  as 


Wsd  ~~  ^sd  /  2)  —  2Zsd  /  d  ,  (20) 

Another  way  to  interpret  the  above  equation  is  that  as  the  disk  is  moved  away  from  the 
observer,  the  corresponding  spatial  frequency  increases  (ie,  narrower  subtense  angle).  From 
Fig,2,  we  see  that  due  to  the  rapid  decay  of  modulation  with  increasing  frequencies,  any  errors 
introduced  by  above  approximation  will  also  decrease  rapidly  when  the  disk  is  moving  away 
from  the  observer. 

Assume  at  this  point  that  the  difference  between  the  target  brightness  and  the  background 
is  small  (but  not  zero),  such  that  (dropping  angular  notations  hereon  for  simplicity) 

X>r~  Lb  or  CL  «  1 ,  (21) 

We  can  relate  the  visibility  contrast  to  the  modulation  transfer  function  via  the  following: 
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(22) 


mwso)= 


^man^ro)  ^min  (V SD  )  j  ,  ^ 

Smx  (WsD )  +  ^min  (y*SD  ) 


=  ( 


4-4w»  _  Q 


4+l, 


-)/Mn  = 


(CL  +  2) 


-/Mn  ~ 


2  Mn 


As  noted  previously  since  both  contrast  terms  depict  the  same  target,  the  informational 
content  is  the  same.  The  difference  exists  only  in  mathematical  forms,  yet  it  gives  a  new 
perspective  about  Secchi  disk  disappearance.  It  is  easy  to  see  from  Gq.  (22)  that  for  C/.=0.02, 
the  above  approximation  introduces  about  1  %  error. 

Rewriting  Eq,  (14), 


D(y)  =  ~ 


ln[tf(y)] 

R 


At  Secchi  depth  Z  =  Zjp,  combining  Eqs.  (19),  (22)  and  (23): 


C-- 


27t8G\ffm 


ln(Jik_). 

2Afn 


(23) 


(24) 


for  the  conditions  0Q  -0.03  [13],  disk  size  t/=0,3m,  for  ZSD  on  the  order  of  10m,  exp(- 
2^^siA>)^  0.002 «L  Let  c=  —  InfC^ /2M0),  then 


C  - - -  — =—  .  (25) 

2jr9Qy/SD  ZSD 

The  above  is  the  derived  horizontal  visibility  range.  For  vertical  Secchi  disk  applications, 
one  can  show  from  Eqs,  (2),  (3)  that 

c  +  K - - - =  —  .  (26) 

Zsd 

Eq.  (26)  is  the  equivalent  form  of  the  basic  equation  for  Secchi  disk  referred  to  hy 
Preisendorfer  [2],  but  now  derived  using  modulation  transfer  theory.  If  we  let 

F  =  £  +  dbf(47T&0 ) ,  then  the  Secchi  disk  depth  (or  visibility  range  in  horizontal  direction) 
can  he  written  as 


cZSD  =  F  for  the  horizontal  case,  and  (27) 

(c  +  K)ZS[>  =  F  for  the  vertical  case.  (28) 

These  have  the  same  form  and  are  directly  comparable  to  the  radiative  transfer  results  |2]. 
Rewriting  Eq,  (27)  gives  the  horizontal  visibility  range  as 

ZSD  =-•  (29) 

c 

Now  let's  consider  a  situation  where  the  black  Secchi  disk  is  used.  Take  the  generally 
accepted  limiting  contrast  as  0.02  from  atmospheric  studies  [  14),  with  M^L  Co=L  c  is  then 
4,6.  Assuming  the  total  scattering  is  0,2  m  1  with  #0=0.03,  then  for  a  30cm  Secchi  disk,  the 
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second  term  in  P  is  given  as  dhtAjrO^c  -0.16 Ic  .  Therefore  the  coupling  constant  F  is 
4.76.  The  horizontal  visibility  range  for  the  black  Secchi  disk  therefore  becomes 

4.8 

visibility - ,  (30) 

c 


which  is  inline  with  previous  studies  [4], 


Fig  3.  Visibility  range  comparisons  between  previous  radiative  transfer  approach  and  current 
MTF  approach.  c(532)  ranges  from  0.24  to  0.43.  Sec  text  for  details. 


Using  absorption  and  attenuation  coefficients  measured  with  an  ac-9  (WetLabs  Inc)  during 
the  GLOW  experiment  (Gauging  Littoral  Optics  for  the  Warfighter,  Sept  17-22,  2001, 
Pensacola,  FL),  calculated  visibility  ranges  using  the  radiative  transfer  and  MTF  approaches 
are  shown  in  Fig.  3.  Visibility  measurements  were  taken  at  lO-foot  (10FSW)  and  30-foot 
(30FSW)  water  depths  (depth  below  the  surface)  each  day,  following  the  same  track 
directions  each  day.  Attenuation  (c)  and  scattering  ( b )  values  at  532nm  are  used  to  represent 
photopic  attenuation  and  scattering  [4]  for  both  approaches.  Each  data  point  represents  the 
averaged  value  at  the  specific  depth  on  each  day*  4.8 ic  is  used  to  derive  visibility  ranges 
representing  previous  studies  [4],  while  Eq.  (29)  is  used  to  obtain  MTF  visibility  ranges  along 
with  typical  parameters  (%=0«Q3,  d=0.3m).  The  two  approaches  generally  agree  well  as 
shown.  However  the  current  approach  gives  a  smaller  predicted  visibility  range  in  dearer 
waters  (c-0.24)  and  a  larger  range  in  more  turbid  waters  (c-0.45).  The  solid  line  depicts  1:1 
ratio. 

3,  Discussion 

To  the  authors'  knowledge,  the  effect  of  disk  size  on  Secchi  depth  has  only  been  implicitly 
embedded  in  previous  theoretical  treatments  of  the  Secchi  disk  method,  usually  with  respect 
to  the  subtense  angle  for  the  eye.  In  other  words,  the  traditional  theory  based  on  the  radiative 
transfer  approach  applies  regardless  of  whether  the  disk  is  1mm  or  10m  in  diameter.  Using 
simple  but  important  assumptions,  this  study  demonstrates  a  similar  Secchi  depth/visibility 
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range  can  be  reached  using  a  more  general  approach  involving  all  spatial  frequencies.  In  doing 
so,  the  Secchi  depth  is  now  explicitly  related  to  the  disk  size  as  well  as  the  scattering.  The 
decay  of  MTF  at  spatial  frequencies  over  the  observing  range  (Fig.  2)  provides  the  straight* 
forward  answer  to  the  Secchi  disk  visibility:  as  the  disk  moves  away  from  the  observer,  the 
spatial  frequency  corresponding  to  the  disk  size  increases,  and  the  MTF  decreases  as  a  result 
of  both  absorption  and  scattering  at  an  increased  rate.  From  Eq.  (19),  we  see  that  the  role  of 
the  total  attenuation  (e)  is  to  dampen  the  MTF  by  the  same  amount  to  all  frequencies,  which  is 
the  basis  of  the  previous  radiative  transfer  approach.  The  current  MTF*based  approach  shows, 
on  the  other  hand,  that  the  scattering  does  not  affect  all  frequencies  equally  Eq.  (19),  (Fig.  2). 
The  only  time  both  approaches  would  converge  is  when  the  second  term  on  the  right  side  of 
Eq.  (19)  becomes  insignificant.  This  is  precisely  the  case  for  the  Secchi  disk  where  the  total 
attenuation  plays  the  dominant  role  at  increased  higher  spatial  frequencies  by  reducing  the 
target  contrast.  This  explains  why  the  original  radiative  transfer  approach  by  Preisendorfer 
and  Dunlley  works  so  well. 

Throughout  the  history  of  Secchi  disk  usage  various  disk  sizes  have  been  used  ranging 
from  a  few  centimeters  to  meters  in  diameter.  The  disk  sizes  most  often  used  are  between  20 
and  40cm  in  diameter;  with  the  30cm  being  the  standard  for  marine  scientists,  while  lake 
researchers  seem  to  prefer  the  20cm  size.  Radiative  transfer  theory  results  do  not  explicitly 
contain  relationships  between  the  Secchi  depth  and  the  disk  size.  Rather,  it  is  embedded 
implicitly  in  the  contrast  threshold  Cv .  From  the  MTF  aspect,  however,  the  functional 
dependence  of  Secchi  disk  visibility  and  disk  size,  and  the  relationship  with  forward  scattering 
Eqs.  (27)  (28),  is  explicit  in  the  formulation.  There  are  other  contributors  to  a  change  in  the 
apparent  size  of  the  disk  such  as  deployment  height  and  viewing  angle,  but  here  we  address 
disk  size  specifically. 

Because  of  the  large  range  of  disk  sizes  used  it  is  important  to  examine  the  impact  disk 
size  has  on  the  overall  visibility  range.  We  can  derive  the  rate  of  variation  of  ZSD  as  a  function 
of  the  disk  size  using  Eq.  (28)  and  demonstrate  that 


dZSp  _  h 

dd  4  M90(c  +  K)' 


(31) 


Applying  typical  values  found  in  coastal  regions  [15}  (ie  a=Q.l4,  ^=0.18  for  bay,  and 
a=0.041,  /?=0.2I  for  coastal  waters,  respectively),  the  change  in  the  range  would  be  3-5%  for 
waters  with  a  small  MSA  value  (0,03),  and  1-1.6%  for  a  larger  MSA  value  (0.1).  Clearly  the 
rate  of  variation  has  a  relatively  small  contribution  to  the  overall  visibility  range,  as  it  only 
affects  a  small  portion  of  t\  Likewise,  using  typical  values  above,  it  can  be  seen  that  the 
scattering  term  represents  roughly  3%  of  the  variation  in  T.  One  may  also  notice  that  the 
approximation  error  of  ^sd  in  Eq.  (20)  does  little  to  the  outcome  of  the  current  theory.  For 
example,  if  a  cycle  of  d  instead  of  dll  is  used,  the  change  in  f  is  only  0.08/c  for  Eq,  (30). 
Together,  these  further  help  explain  the  convergence  of  radiative  transfer  and  the  MTF 
approaches. 

Blackwell  also  showed  |6|  that  the  limiting  contrast  Q  is  a  function  of  disk  size  (or 
angular  subtense)  and  adaptation  luminance,  i v.dgfdd  ^0,  Analysis  with  the  radiative 
transfer  method  has  shown  a  small  (<i%)  change  in  range  visibility  due  to  disk  size  changes 
(43cm  to  60cm)  1 1  ].  This  effect  can  also  he  explained  by  Fig. 2.  From  the  modulation  transfer 
perspective  it  can  be  seen  that  the  contrast  decay  of  the  spatial  frequency  reaches  a  plateau  at 
higher  spatial  frequencies.  In  other  words,  due  to  the  fast  decay  in  the  MTF  it  does  not  really 
matter  much  if  the  disk  size  is  20cm  or  30cm,  This  is  an  impact  only  when  Z$&  is  very  small 
or  when  the  disk  is  very  large  (ie  low  frequencies).  This  rapid  decay  is  the  key  reason  behind 
observations  that  the  target  size  does  not  significantly  alter  the  disappearance  range  of  the 
disk.  The  further  away  the  disk  is,  the  higher  the  spatial  frequencies  it  corresponds  to,  and 
where  the  MTF  at  these  frequencies  flattened  out  is  about  the  same  for  the  30cm  and  the  20cm 
disk. 
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As  for  disk  color*  from  the  modulation  transfer  theory  Eqs,  (II)  (25),  and  everything  else 
being  equal,  the  usage  of  the  black  disk  versus  black  and  white  alternating  disk  should  have 
the  same  outcome  since  they  both  have  M<pl,  This  has  been  confirmed  by  daily  Secehi  depth 
measurements  over  a  6  month  time  span  done  in  Skagit  River  [8].  Therefore,  both  types  of 
disks  are  equally  useful  when  measuring  underwater  visibility  and  water  quality  [4], 

The  contrast  term,  is  somewhat  less  defined  since  it  is  solely  dependent  upon  each 
observer’s  ability  to  discern  contrast*  and  subsequently  is  also  a  function  of  conditions  at 
which  measurements  are  taken  (ie,  ambient  light*  sea  state,  surface  glint,  sun  angle,  shadow  or 
sunny  side  of  vessel,  disk  types*  bottom  types,  and  so  on).  However*  claims  have  been  made 
to  indicate  very  small  variations  (<l%)  among  different  observer  groups,  volunteers  versus 
professionals  amongst  lake  monitors  in  Florida  [16],  This  result  seems  inconsistent  with  the 
theory  which  suggests  that  for  differing  background  adaptation  luminance,  combined  with 
different  contrast  values  (0,02  to  0.005),  the  possible  f  values  can  range  from  5,9  to  9.3  for  a 
white  disk  with  85%  reflectance*  even  under  optimal  observation  conditions*  per  Preisendorter 
[2],  Using  typical  coastal  region  parameters  quoted  earlier*  the  results  here  show  V  ranging 
from  6.5  to  10. 1  Eq.  (27)  under  similar  conditions  as  Preisendorter  [2J,  Other  studies  report 
different  values,  such  as  r=8.68  [17],  and  r=4,8  [4],  As  expected  the  V  value  can  be 
significantly  higher  under  conditions  of  strong  scattering.  For  example,  when  h- 2  and  all 
other  parameters  stay  the  same  as  in  the  earlier  black  disk  case  (where  T =4.76  (6=0,2))  then 
T=6,2  will  be  obtained  with  the  current  formulation.  The  larger  range  of  f  values  can  be 
explained,  in  pari,  by  the  differences  in  the  scattering  contribution  to  the  total  attenuation. 

The  modulation  transfer  method  presented  here  has  the  benefit  of  a  general  approach,  and 
is  applicable  to  other  underwater  visibility  issues  including  self-illuminating  targets.  The 
disappearing  frequency  corresponds  to  any  contrast  related  features,  such  as  patterns  and 
textures*  so  long  as  the  dependant  spatial  frequency  is  not  so  small  that  it  is  limited  by  the 
system  hardware  before  the  medium  Eq.  (9),  The  only  difference  when  applying  the  theory  to 
different  targets  is  in  the  threshold  contrast  that  is  required. 

The  small  angle  approximation  by  Wells  is  applicable  to  most  vision  and  imaging  systems 
with  a  field-of-view  cone  less  than  20°,  It  considers  photons  scattered  within  this  angle  as  pan 
of  the  signal  to  be  included.  This  inherently  assumes  strong  forward  scattering  in  the  medium. 
Monte  Carlo  simulation  results  indicate  that  at  the  typical  spatial  frequency  for  a  Secchi  disk 
range  (>50  cycles  per  radian),  the  small  angle  approximation  error  is  about  2%  for  up  to  10 
attenuation  lengths  [15],  in  addition*  studies  also  show  that  the  exact  shape  of  the  volume 
scattering  function  Eq.  (18)  does  not  affect  the  outcome  significantly  {<!%)*  as  it  has  little 
impact  on  the  cumulative  effect  of  PSF*  and  therefore  the  MTF  [13].  However,  departure 
from  the  small  angle  approximation  assumption  will  invalidate  the  modulation  transfer  theory 
applied  here.  For  instance,  the  theory  is  not  valid  for  a  Raleigh-type  scattering  medium. 

The  assumption  to  eliminate  the  exp(-2^iA)  term  in  Eq.  (24)  may  be  invalid  in  very 
turbid  environments.  For  example,  if  visibility  is  only  1.5m,  ^sd-  1 0  and  exp(-2jn|fsn(y  -0.15* 
in  which  case  the  contribution  should  he  included  as  it  represents  -  15%  variation.  However, 
from  discussions  above,  since  this  entire  term  contributes  to  only  -3%  to  f  overall,  the 
conclusions  arrived  at  in  Eqs.  (26) -{29)  still  hold  well.  Further,  a  closer  examination  at 
different  spatial  frequencies  reveals  that  this  term  would  affect  only  low  frequencies  (Fig.  2), 
and  therefore  is  not  a  major  concern  in  the  case  of  the  Secchi  disk. 

Modulation  transfer  decay  (blur)  due  to  the  turbulence  of  the  medium  is  not  accounted  for 
in  our  discussion,  nor  is  other  factors  that  reduce  visibility  such  as  surface  glint  and  capillary 
waves.  In  addition,  one  cannot  ignore  an  implicit  assumption  when  applying  modulation 
transfer  theories,  namely  the  shift  invariance  of  the  system.  The  shift  invariance  roughly 
translates  into  requirements  in  the  field  that  Secchi  disk  inradiance  does  not  change  over  the 
period  of  observation.  Fortunately,  the  human  vision  system  automatically  adapts  and 
eliminates  many  of  these  effects.  This  requirement  is  important,  however*  when  electro-optic 
imaging  systems  are  used  or  designed  based  on  the  above  discussed  theory. 
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4.  Conclusion 


The  disappearance  of  the  Secchi  disk  is  the  result  of  the  reduced  modulation  contrast  of  the 
disk  in  the  water  medium.  This  is  brought  about  by  absorption  and  scattering,  with  the  later 
playing  a  more  significant  role  due  to  its  effect  on  the  contrast  reduction.  We  demonstrate  that 
by  applying  the  modulation  transfer  theory,  the  Secchi  depth  and  what  is  generally  referred  to 
as  the  horizontal  visibility  can  be  found  using  a  more  general  approach  involving  all  spatial 
frequencies.  We  show  that,  when  typical  optical  properties  are  used,  predicted  visibility 
ranges  are  comparable  to  those  obtained  using  radiative  transfer  methods.  This  result  is 
obtained  under  the  assumption  that  the  small  angle  scattering  approximation  is  valid,  which 
demands  strong  forward  scattering  as  the  pre-eon dition.  The  effect  of  disk  size  to  Secchi  depth 
is  included  in  the  current  form,  both  implicitly  embedded  in  contrast  threshold  as  previous 
studies,  but  also  explicitly  expressed  as  well.  The  improvement  in  this  theoretical  work  is  that 
not  all  spatial  frequencies  (ie  sizes)  are  equally  critical  in  visibility  studies.  In  the  special  case 
of  the  Secchi  disk  it  is  shown  that  the  modulation  transfer  approach  and  the  radiative  transfer 
approach  converge.  This  is  due  to  a  rapid  decrease  in  the  MTF  and  a  flattening  in  the  high 
frequency  region  due  to  attenuation  of  transmitted  contrast.  The  MTF  theory  applied  in  this 
work  provides  a  different  perspective  from  which  a  new  understanding  about  underwater 
visibility  issues  can  be  achieved,  and  from  which  new  applications  of  Secchi  disk  visibility  to 
target  identification  and  field  instrumentation  can  be  envisioned. 
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